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Abstract. We study the nonequilibrium properties of an electronic circuit composed 
of a double quantum dot (DQD) channel coupled to a quantum point contact (QPC) 
within the framework of stochastic thermodynamics. We show that the transition rates 
describing the dynamics satisfy a nontrivial local detailed balance (LDB) and that the 
statistics of energy and particle currents across both channels obeys a fluctuation 
theorem (FT). We analyze two regimes where the device operates as a thermodynamic 
machine and study its output power and efficiency fluctuations. We show that the 
electrons tunneling through the QPC without interacting with the DQD have a strong 
effect on the device efficiency. 


1. Introduction 

Semiconducting multichannel circuits made of quantum dots and quantum point 
contacts are nowadays commonly devised and studied experimentally [1, 2, 3, 4, 
5, 6, 7, 8, 9, 10]. The progress in the control of electronic temperatures at the 
meso-scale [7, 8] has for instance driven the experimental [11, 12] and theoretical 
[13, 14, 15, 16, 17, 18, 19, 20, 21, 22] study of their thermoelectric properties. In the 
isothermal case, these circuits have also been used to probe the fluctuating properties of 
heat and matter transfers using counting statistics experiments [5, 23, 9, 10]. A circuit 
of particular interest in that regard is the DQD channel probed by a QPC detector. 
It has been used to perform the bidirectional counting of single electrons in the DQD 
channel [24, 23] and to provide the first experimental verification of the current FT in 
mesoscopic physics [25]. Theoretically, several studies have analyzed the backaction of 
the QPC detector on the mean current of the DQD channel, as well as on its statistics 
[26, 27, 28, 29]. The QPC detector was also shown to modify the thermodynamic affinity 
of the DQD channel while preserving the FT symmetry in the DQD circuit [27]. In this 
latter work, the tunneling events in the QPC were treated non perturbatively to account 
for possible high transparency and as a result, the combined DQD-QPC statistics was 
not accessible within this approach. 

In this paper, we study the nonequilibrium thermodynamics of the DQD-QPC 
circuit using stochastic thermodynamics [30, 31, 32]. We consider the general case 
where the QPC and DQD reservoirs may be at a different temperatures and chemical 
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potentials. The counting statistics of the energy and matter currents across both 
channels is calculated using the modified quantum master equation formalism [33] 
assuming weak coupling between the DQD and its reservoirs as well as between the 
reservoirs composing the QPC. Even though in the isothermal case this last assumption 
may seem more restrictive compared to the nonperturbative approach of Ref. [27, 34], 
it enables us to analytically calculate the join distribution of the energy and matter 
currents across both channels, and to identify the entropy flows associated to the 
exchange processes at hand. As a result, we are able to derive a bivariate FT for 
the statistics of the energy and matter currents in both channels. 

An interesting feature of this setup is that the microscopic processes associated 
to transitions in the DQD involve more than one reservoir at a time. In particular, 
transitions in the DQD induced by the QPC have corresponding transition rates 
proportional to the product of Fermi functions in both reservoirs of the QPC. The 
microscopic processes underlying such transitions involve the tunneling of an electron 
between the QPC reservoirs which exchanges a fraction of its energy with the DQD. As 
a result, the transition rates cannot be written anymore as a sum where each term only 
involves one single reservoir. Despite this non additivity of the rates, the local detailed 
balance (LDB) [33, 35] is shown to hold, and is explicitly written in terms of the fluxes 
of entropy from the reservoirs involved in the transitions. 

Furthermore, since some electrons tunnel between the QPC reservoirs without 
exchanging energy with the DQD, they do not induce transitions in the DQD but need 
to be taken into account in the counting statistics. We show that their statistics is well 
described by the modified quantum master equation and corresponds to the Levitov- 
Lesovik formula [36, 37, 38, 39, 40] to second order perturbation theory in the coupling 
between the QPC reservoirs. Though such processes do not reveal themselves in the 
master equation for the DQD populations, they are shown to strongly impact the circuit 
performance. 

We proceed by analyzing two different regimes where the circuit operates as a 
thermoelectric and as a current converter, respectively. We identify the optimal working 
conditions to reach large average output power and high macroscopic efficiency and study 
the statistical properties of the output power as well as of the efiiciency, as recently 
proposed in Refs. [41, 42, 43]. 

The paper is organized as follows: The model is introduced in section 2. In section 
3, we derive the modified quantum master equation to calculate the counting statistics 
of energy and matter currents, and analyze the microscopic processes contributing to 
the transition rates. In section 4, using the LDB property of the rates, we identify the 
entropy flows associated to each microscopic processes in the circuit, and the steady 
state FT is derived. The expressions for the average energy and matter currents as well 
as the average irreversible entropy production are also provided. The thermodynamic 
analysis of the circuit operating as a thermoelectric and as a current converter is done 
in section 5. Conclusions are drawn in section 6. 
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The DQD channel is made of two quantum dots A and B, each connected to its own 
reservoir, labeled by j = 1 and 2 respectively. The QPC is the junction between 
reservoirs j = 3 and 4. The circuit is drawn in Fig. 1. 



Figure 1. Double channel circuit made of a DQD and a QPC. Electrons can be 
exchanged between reservoirs 1 and 2 across the DQD, and between reservoirs 3 and 4 
of the QPC. There is no electron transfer between the two channels, but the electrons 
tunneling through the QPC are affected by the charge state of the quantum dots via 
Coulomb interaction (wavy lines). 

The Hamiltonian of the circuit is given by 
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j = \ 

(1) 

where i^oQD denotes the DQD Hamiltonian, 

(2) 


k 


is the reservoir j Hamiltonian expressed in terms of the creation (anihilation) operators 
^jk of the reservoir single-particle states with energy and V is the interaction 
Hamiltonian between the DQD and the reservoirs. 

The DQD Hamiltonian is given by 

-f^DQD = GaCaCa + escj^CB + T [c\cb + C^jCa) = ns|s)(s|, (3) 

expressed in terms of the single-dot annihilation (creation) operators in each QD, ca/b 
(4/s)’ single dot states with energies e^/^, and of the tunneling amplitude T between 
the two dots. This Hamiltonian can alternatively be expressed in terms of the many- 
body eigenstates |s) of i^DQD with energies Eg. In the following, we consider a regime 
in which the sum in the last term of (3) can be restricted to the empty eigenstate |0) 
and the single-occupied eigenstates |-|-) and |—) (see the Appendix). This assumption 
is often justified at low temperature since multiple charging in one of the quantum dots 
requires large amounts of energy. 
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The interaction Hamiltonian can be split into 
2 

^ = + (4) 

J = 1 

The first term is the sum of the tunneling Hamiltonians between the DQD and the 
reservoirs 

^- = XI Z] {c\cjk + cJfcCd) = X X 

d=A,B k s=-\-,— k 

with the tunneling amplitudes and in the DQD single and many-body basis, 
respectively, and for j = 1 and 2. Each dot is only connected to its own reservoirs, 
i.e. t^Q = t 2 yi = 0- This part of the Hamiltonian is responsible for the charging and 
discharging of the DQD through the exchange of electrons with reservoir j = 1 and 2 . 
The QPC reservoirs 3 and 4 are directly coupled through the tunneling Hamiltonian 

V'34 = E^“'(4-:4f+4.C3t). ( 6 ) 

kk' 

where denote the bare tunneling amplitudes between the QPC reservoirs, i.e. 
independently of the DQD state. 

Finally, electron tunneling in the QPC is also sensible to the charge state of the 
DQD due to the Coulomb interaction. This interaction is modelled by the capacitive 
coupling 

U4^^ = X {^a'caCA + (4^4^' + 4fc'C3fc) (7) 

kk' 

= EE( (c\^c^k' + 4fc'C3fc) (8) 

kk' 

in terms of the capacitive couplings t^A/B tunneling amplitudes . Expression 

( 8 ) shows that some electrons in the QPC can induce transitions in the DQD, exchanging 
energy with the DQD while tunneling between reservoirs 3 and 4. 

The Hamiltonian term describing the bare tunneling in the QPC, V 34 , is included 
in the interaction V which will be subsequently treated to second order in perturbation 
theory. This is in contrast to our previous work [27] where the bare tunneling in the QPC 
was treated non-perturbatively. The present approach has the advantage to treat all 
the energy and matter transfers on the same footing, allowing us to develop a consistent 
thermodynamic description of the full circuit and to evaluate the joint full counting 
statistics (ECS) of the energy and matter currents in both DQD and QPC channels. 

3. Counting statistics 

The ECS of the fluxes in weakly coupled open quantum systems can be calculated using 
the modified quantum master equation formalism [33] . The statistical properties of the 
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energy and matter cnrrents flowing ont of the reservoirs and integrated over a time t, 
/S.Ej and AA/^ , are determined by the generating function (GF) 

Xj, t) = 


The probability distribution of the fluctuating energy and matter fluxes = —AEj/t 
and = —ANj/t is then obtained by applying an inverse Fourier transform to the GF 

(9) 

/*27r 




n* 


dJl 

271 


L J 




271 


L J 


_l\.^t). ( 10 ) 


In the following, we evaluate the GF (9) by performing the FGS of the energy 
and particle number operators within each reservoir, respectively given by ( 2 ) and 
Nj = Following [33], we introduce the modified Hamiltonian 

Xj ) = e * ^^"=1 ^ He~^ ^^9 ( 11 ) 


where the counting parameters and Aj for j = 1,..., 4 keep track of, respectively, the 
energy and matter fluctuations in the reservoirs. The DQD and reservoir Hamiltonians, 
-f^DQD and Hj respectively, remain unchanged after this transformation. However, the 
interaction Hamiltonians transform according to 

|s> (0| Cj, + h.c.) , (12) 

5 =+,— k 


A,) = + h.c.) , (13) 

kk' 

and 

r,?«A{.A,) = EE (^r,^/|s)(s'|) + h.c.) . (14) 

ss' kk' 

With these definitions, the GF can be expressed as 

Gi^j, Aj, t) = Tr t)} , (15) 

where the modified density matrix p(^j,Xj,t) satisfies the modified quantum master 
equation 

idtpi^j, Xj, t) = H{^j/2, Xj/2)p{ij, Xj, t)-p{Cj, Xj,t)H{-^j/2, -Aj/2).(16) 

The initial density matrix p(0) of the total system is assumed of the factorized 
form p(0) = ps(0) Hj ®Pj, where ps(0) is an arbitrary DQD density operator and 
Pj = exp {—l5j{Hj — PjNj — (f)j)} denotes the grand-canonical density operator in the 
reservoir j with inverse temperature (5j = (k-BTj)~^ and chemical potential pj with 
j = 1,...,4. The corresponding thermodynamic grand-potential is denoted <pj = 
In [Tr {exp — PjNj)}}]. This factorization assumption has no implication 

because only steady state properties will be considered in the following. 

In the weak coupling limit to the reservoirs, where the interaction parameters 
Tjg, and are assumed small enough, the dynamical equation (16) leads to a 
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closed modified quantum master equation for the reduced density matrix of the system 
[44, 45, 46, 47, 48] 

Psi^jAjA) = TrK{pi^j,Xj,t)} , (17) 

where TrR{-} denotes a trace over the reservoirs Hilbert space. ps{0,0,t) is the DQD 
reduced density matrix. 

A common assumption in the present context is that the DQD free oscillations, 
characterized by the frequencies cUgs' = Eg — Eg', are fast compared to the relaxation 
time scale tr induced by the reservoirs on the DQD. One can then apply the rotating 
wave approximation (RWA) [49, 50, 51, 52, 53] which consists in an average of the system 
free oscillations over a time scale At which is intermediate between 

Tc -C At -C tr, (18) 

where tc denotes the short correlation time in the reservoirs. As a result, the effective 
dynamics of the DQD populations, ^^(Q, Aj, t) = {.s\ps{iijA^ji ^)|'S), and the coherences, 
{s\ps{iij, i^j,t)W) for 5 7 ^ s', decouple. 



{(iA,} = {±l,0,0,0} 


{5iVQ = {0, ±1,0,0} 


a) 


{5Ej} = {±a;so, 0,0,0} 


b) 


{<5AQ = {0,±a;,o,0,0} 



{SNj}^ {0,0, ±1,^1} 
c) {SEj}^{0,0,±e,T{e-uj+.)} 


QAQ = {0,0,±1,t1} 

d) {(iAQ = {0,0,±(e-a;+_),Te} 


Figure 2. Illustration of the several microscopic processes inducing transitions in the 
DQD. Each sub-figure illustrates pairs of processes which are time-reversed of each 
other. The vectors {SEj} and denote, respectively, the energy and particle 

number changes in the reservoirs associated to each microscopic processes. 

Under the aforementioned hypotheses, the diagonal elements of the DQD reduced 
density matrix satisfy a Markovian master equation of the form 

dtSiijAj.t) = W(5,,A_i) ■g(5,,A,.t), 


(19) 
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where we introduced the vector notation 


^3 A) 




( 20 ) 


\ 9- {Cj ; ) t) J 

with the matrix product denoted by and where the counting parameters dependent 
rate matrix W(^j, Aj) is expressed here below in terms of the transition rates between 
DQD states. 

The reservoirs j = 1 and 2 induce random charging and discharging of the DQD due 
to the tunneling interaction (5). The corresponding microscopic processes are depicted 
on Figs. 2 a) and b). During such tunneling events, the particle number in reservoir 
j = 1 or 2 changes by an amount 5Nj = 1 {SNj = — 1) when it charges (discharges) 
the DQD. On the other hand, the energy change in the reservoir can take the values 
5Ej = ±C(Jso, depending on which many-body state |s), with s = + or —, is involved in 
the transition. The charging and discharging rates induced by reservoirs j = 1 and 2 
are given by 


Ojg Y'jgfjiiOgoj and bjg rj\ 5 (l fj(ujgQ^^ (21) 

for s = + or —, in terms of the Fermi distribution of single particle states in the reservoir 
j, fj{x) = (1 + exp I5j{x — and of the rate constants 

27r 27r 

( 22 ) 

k 

We took the continuum approximation for the electron density of states in the reservoirs, 
denoting the energy-resolved tunneling amplitudes by Tjs{e), and the density of electron 
states by Pj{e). 

On the other hand, the QPC also induces transitions between the DQD states. 
Though there is no exchange of electrons between the QPC and the DQD, electrons 
tunneling between the reservoirs j = 3 and 4 may exchange energy with the DQD 
(mainly through photon exchange [6]), thus driving transitions between the single- 
charged states |-|-) and |—). These processes are illustrated in Figs. 2 c) and d). The 
corresponding transition rates are given by 


Oi'(e) = riy(e) /i(e) (1 “ /i'(e “ ‘^+-)) 

= rii'(e) (1 - /i(e)) fj'{e - uj+_) (23) 

for jj' = 34 and 43 where the energy-dependent rate constants Tjji (e) are given by 

/j /jr 

r34(e) = - esk)Sie - a;+_ - (24) 

kk' 

Ajr 

= |r-+(e, e - w+_)|V3(e)p4(e - ^^+-) (25) 

/j /jY 

r43(e) = - ^+- - Ofc)<^(e - OfcO (26) 

kk' 

Zj-TT 

= ■j^|T+-(E-‘^+-.E)l"P3(E-i^+-)P4(£) (27) 

( 28 ) 
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in terms of the energy resolved tunneling amplitudes Tag'^e^e'). Interestingly, the 
transition rates (23) are written as a product of Fermi functions in both QPC reservoirs 
and as such, cannot be written as a sum of individual reservoir contributions. 

The tunneling events between the DQD and the reservoirs j = 1 and 2 contribute 
to the rate matrix through the matrix elements 

|Wte,A,)],„= (29) 

i=i,2 

[W(?„A,.)]„. = Y. (30) 

i=i.2 

where the counting parameters and Xj keep track of the net fluxes of energy icUgo 
and particles ±1 flowing out of reservoir j at each such transition. 

The contribution from the QPC transfers can be separated into two categories, 
depending on whether or not the electrons tunneling between reservoirs 3 and 4 exchange 
energy with the DQD. In the first case, tunneling events in the QPC contribute to the 
rate matrix through the components 

1W(«J.A,)]^._ = Y.( *c„-(e)efe-«>''™+-)+^J-V) (31) 

jj' 

[W(^J,A,)]_^ *(i,,-(e)e(-«J«+«/<-“+-l-^i+V) (32) 

jj' 

where the sum in these last two equalities runs over jf = 34 and 43. Each transition 
involves a net transfer of one electron from reservoir j to reservoir / of the QPC or vice 
versa. If the energy of the outgoing electron coming from j is e, it enters reservoir j' 
with energy e ± (^+_ depending on whether it emits or absorbs energy from the DQD. 


I-) 1+) 

• • 



{(iEQ = {0,0,±e,Te} 

{5iV,} = {0,0, ±1,^1} 


Figure 3. Illustration of the tunneling processes in the QPC that do note induce 
transitions in the DQD. 

Finally, the contributions from the electrons tunneling in the QPC without 
exchanging energy with the DQD appear along the diagonal elements of the rate matrix 

1W(6, A,)]„ = - ^ [W(0,0)].,, + G.iii, Aj). (33) 

s'^s 
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The first term in the right-hand side of this equation ensures the conservation of the 
probability for the occupation probabilities in the DQD when the counting parameters 
are set to zero while the second one accounts for the tunneling of electrons through the 
QPC without interaction with the DQD (see Fig. 3). As a matter of fact, Gs{ij,Xj) is 
the GF of the energy and particle transfer in the QPC given that the DQD is in state 
|s), i.e. 

Gs{ijAj) = j de'jsie) 

[/ 3 (e)(l - /4(e)) (1 - e(^3-A4)gde3-e4)) +/,(e)(i _ f^(e)) (l - . (34) 

It turns out that this is the Levitov-Lesovik formula [36, 37, 38, 39, 40] to second order 
in the tunneling amplitude Ts{e), or equivalently, to first order in 

7 s(e) = 27r|I/3s(e)|V3(e)p4(e), (35) 


which results from the fact that we treated the interaction (8) perturbatively. In 
contrast to many previous work on transfer through quantum dots (such as [54]), where 
each microscopic process is associated to an actual transition in the QD, the present 
circuit provides a nice example of an open quantum system in which the microscopic 
processes do not necessarily affect the system (DQD) populations while contributing 
to the energy and particle flows out of the reservoirs. The modified quantum master 
equation formalism is thus essential in order to keep track of these processes. 

The Fourier transform of the diagonal matrix elements (20) 


p{AEj, ANj,t) = 


n 


271 


p27r 


jj- dXj 


27T 


L J 




gives the join probability distribution [p(A£), ANj,t)]g = Ps{AEj, ANj,t) of observing 
the system in state |s) at time t and the changes in energies AEj and particle numbers 


ANj in each reservoir. 


By applying a Fourier transform to the modified rate equation (19), we get a master 
equation describing the dynamics of the DQD as well as the exchange processes with 
the reservoirs 

dtp{AEj, ANj, t)= [Y[ [ddEj] [dSNj] W{SEj, SNj) ■ p{AEj - 5Ej,ANj - SNj,t), (37) 


where the rate matrix W((j£'j, dNj) is obtained as the Fourier transform of the modified 
rate matrix 


W(SE,J/Vj) = 


n 


L J 


277 


r27r 


jj- dXj 


271 


w(-iQ, -iXj). (38) 


The Fourier transform of the modified rate matrix is easily taken by using the relations 


r>27r 


277 


= d 


a,0 


and 


/ 


dx . 


27r 


= h(«) 


(39) 
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in terms of the Kronecker delta symbol and the Dirac delta distribution (5 (q;). 
Accordingly, the rate matrix 'W{5Ej,5Nj) is obtained by making the following 
substitutions 


e=*=b"j S{SEj T ctj) and e^^^' ^ (40) 

in the modified rate matrix elements (29)-(33). 

By integrating equation (37) over the energy and particle fluctuations AEj and 
ANj, or equivalently by setting the counting parameters to zero in the modified rate 
equation (19), we obtain a stochastic master equation for the occupation probabilities 
in the DQD 


p(i) = W . p(t) 


(41) 


with rate matrix given by 

/ —Oi-i- — 0(2-1- — — ®2— ^1-1- E ^2+ 

ai-(- -|- 02+ ~bl+ — b2+ — 0^34 ~ ^43 

tti— 02— d$4 C?43 

and where 


W 




bi— 62 — 

C 34 -|- C 43 

-C34 - C43 - bi - - 62- J 


,(42) 




= J decjpie) 


and 


d„, - j 


(43) 


Now, by formally solving the rate equation (19), the GF of the energy and matter 
currents can be expressed as 

Xj, t) = l- (44) 


where 1 denotes the line vector (1, 1 , 1) while po is the initial occupation probability of 
the DQD states. All the moments of the currents can be obtained by taking multiple 
derivatives of the GF with respect to the counting parameters. 

At steady state, the statistics of the currents is captured by the cumulant generating 
function (CGF) 

, ^i) = - lim 7 In , Xj , t). (45) 

t—)-oo t 

Using expression (44) for the GF, one sees that the CGF is obtained as the dominant 
eigenvalue of the rate matrix W(^j,Aj), independently of the initial condition on the 
system po. 


4. Nonequilibrium thermodynamics 

4 . 1 . Local detailed balance and fluctuation theorem 

The charging and discharging rates (21) depicted on Figs. 2 a) and b) satisfy the LDB 
condition [35] 
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in terms of the inverse temperature 13j and chemical potential fij of the reservoir involved 
in the transition. This property is a direct consequence of the Kubo-Martin-Schwinger 
(KMS) condition which is satished by the equilibrium correlation functions of the 
reservoirs [33]. 

Regarding the QPC induced transitions on the DQD, it has previously been noted 
[27, 34] that the corresponding total rates do not satisfy a LDB, unless the QPC is 
assumed to be at equilibrium {(3^ = /54 and = 114 ). However, we showed in the 
previous section that in the weak coupling limit, one can identify the pairs of microscopic 
processes related by time reversal (cf. Figs. 2 c) and d)), and write the total rates as 
a sum of contributions from such elementary processes, see equations (31) and (32) 
together with (23). Each of these contributions satishes the LDB condition 

- Pi) + ^ 3 ' (e - - Pi') > (47) 

«ii'(ej 

where Cjj/ (e) is the rate at which electrons with energy e tunnel from reservoir j to j' in 
the QPC while releasing an amount of energy to the DQD, and djp{e) is the rate 
of the associated time-reversed process. 

We note that the right-hand sides of expressions (46) and (47) is nothing but the 
entropy flowing from the reservoirs during these processes 

Ass'S{SEj,SNj) = -J2 (^®) 

j 

and where the components of the vectors 


{SEj} = {SEi, SE 2 , SEs, 6 E^} and { 6 Nj} = SN 2 , SNs, SN^} 


(49) 


denote the changes in energy and particle number in the reservoirs associated to each 
microscopic transition, as given in Fig. 2. Relations (46) and (47) can then both be 
rewritten in terms of the rate matrix (38) as 


In 


W(5E„SN,) 


W{-SEj, -SNj 


= A„,S(SEj,SNj), 


(50) 


where the transition rate W (—5Ej, —SNj) 

of the one associated to 'W{5Ej,5Nj) 
transition rates (29)-(33) satisfy 


corresponds to the time-reversed process 
Relation (50) implies that the modified 


1W(0, V)]„, = [W(ft - ij, - A,)],,,. (61) 

These relations also hold for the diagonal elements of the rate matrix, given in equation 
(33), Gsi^j, Xj) = Gs{f3j - ij, -PjPj — Xj). 

They ensure the invariance of the characteristic polynomial of the matrix W , Xj ) 
under the transformations 


^3 1^3 ~ ^3 


and 


Xj -/3jttj - Xj. 


(52) 
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Since the CGF is obtained as the largest eigenvalue of the modified rate matrix 
W(^j, Aj), this invariance property leads to the symmetry 

Q{ij, Aj) = g{l5j - -pjfij - Xj). (53) 

We further note that the CGF only depends on the differences 

6 - ^ 4 , 6 “ ^ 4 , 6 - ^4 and Ai - A 2 , A 3 - A 4 . (54) 

This can be directly shown by verifying that the characteristic polynomial of the rate 
matrix W(^j, Xj) only depends on the differences (54). It is a direct consequence of the 
conservation of the total energy and particle number. In addition, the particle numbers 
in each sub-channel is also conserved since electrons cannot tunnel between the DQD 


and the QPC. Accordingly, we introduce the CGF 
^(6)-^!)-^ 3 ) = "^i)l^4=A2=A4=0 • 

where I = 1, 2, 3. The CGF (55) satisfies the steady state current FT symmetry [55, 56] 
Qi^i, Ai, A3) = g{A^E - 6, Al,-X^,A%- As) (56) 

in terms of the thermodynamic forces applied to the system: 

A^^ = 134 — I3i for / = 1, 2, 3 (57) 

A\f = PiUi - 132112 and = /^s/xs - / 54 /r 4 - (58) 


More explicitly, this FT can be expressed in terms of the join probability 
distribution 

p(4,4,4,t) = / ^4 y dj% J d4p(4,4,t) (59) 

as 


lim - In -^(4; 4^ 4 A) 
t^oot 


3 


E A^ 7^ -l qi 71 -L 43 t 3 
l=l 


(60) 


which makes explicit reference to the steady-state entropy production (right-hand side) 
generated by the fluxes against the thermodynamic affinities (57)-(58). 

In the isothermal setups, one obtains a bivariate FT for the join distribution of 
particle currents through each channel 

lim ^ In f ^4 “ /^ 2)4 + / 5(/^3 - 11a)J%, (61) 

thus extending previous result obtained in [27] , where the statistics of the current in the 
QPG was not assessed. 
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Using (45) together with (44), we can formally write the steady-state mean energy and 
particle currents out of reservoir j as 

(Ji,) = - lim = 1 • %W(0,0) • P (62) 

t—)-oo t 

(4) = - lim = 1 • 4W(0,0) ■ P, (63) 

t^OO t 

where P = hmf_j,oo po denotes the vector of steady-state occupation probabilities 
Pg = [P]^ on the DQD. The steady-state probabilities are directly obtained by solving 
the equation W- P = 0. 

Using the rate matrix of the present model W given by (42), we find for outgoing 
currents of particles and energy from reservoirs j = 1 and 2 that 

(4) = XI 4*^0 - bjgPg) , (4) = X {ajsPo - bjgPs). (64) 

S = + ,- 5 = + ,- 

The particle currents for reservoirs 3 and 4 can be expressed as the sum of two 
contributions 

{J%)u = J de{JN{e))u, for u = d,i, (65) 

in terms of the energy resolved currents 

(^v(e))i = ^•^^(e) - / 4 (e)) ( 66 ) 

{JN{e))d = (c34(e) - C 43 (e))P_ -F ( 0 ^ 43 ( 6 ) - d 34 (e))P+. (67) 

The current (Jat)* is conveyed by electrons which tunnel between the QPC reservoirs 
without affecting the DQD while (Jjv)d is the current of electrons that induce transitions 
between states |-|-) ^ |—) in the DQD through the exchange of an amount ±a;+_ of 
energy with the DQD. The energy currents out of the reservoirs 3 and 4 are in turn 
given by 

( 4 ) = j dee {{.JN{e))d + (Jiv(e))*) + <^+- (C 43 P- - (^434-) (68) 

(4) = - J dee ((Jiv(e))d+ {JNie))i) + lo+- {cs^P- - c?34P+) • (69) 

From these expressions, one readily verifies the conservation laws 

^^(4) ~ 0) (4) = “(4) ^^6 (4) = {jff)d + (4)i = “(4)- ('^ 0 ) 

i=i 

The entropy flow from the environment reads 

= ( 71 ) 

i=i 

while the rate of entropy production Si in the whole setup can be expressed as 

S. = d,S+{Js.), 


(72) 
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where S = —'^^Pslnps denotes the Shannon entropy of the system whose time 
derivative vanishes at steady state. As a result, the steady-state rate of entropy 
production can be linked to the energy and matter currents by 

4 

s, = (*.) = - M)). (73) 

i=i 

Using the conservation laws (70), this equation can be rewritten as a sum of 
terms that can be interpreted as the dissipation generated by each current against its 
thermodynamic affinity 

3 

Si = ^e{Je) + ^n{Sn) + 

i=i 

This expression takes the same form as the left-hand side of equation (60), which can 
be used to prove that it is always non-negative, i.e. Si > 0. The entropy production 
(74) plays a central role for the thermodynamic analysis of our device driven out of 
equilibrium by thermal and chemical potential gradients. It is essential in order to 
define proper notions of efficiency when the device is set to work as a thermodynamic 
machine. In the next section, we make use of our analysis in order to characterize the 
performance of our device when operating as a thermoelectric and an isothermal electric 
converter. 

5. Device operating as a thermodynamic machine 

5.1. Regime of thermoelectric conversion 

We consider a regime in which the QPC is the hot reservoir with inverse temperature 
Ph = ^3 = Pa whereas the DQD reservoirs j = 1 and 2 constitute the cold reservoir with 
inverse temperature Pc = Pi = P 2 - A fraction of the heat flow from the QPC 

Q = J% -\- (75) 

may then be converted by the DQD into electro-chemical work against a bias Afi = 
III — p 2 applied between the reservoirs 1 and 2 

TV = -A/r4. (76) 

We assume a vanishing bias in the QPC so that the working regime of our heat engine 
is 

Ph < Pc, R3 = Ra and A/r < 0. (77) 

The irreversible entropy production (74) reduces to 

= -Pcm + {Pc - Ph) {Q) > 0, (78) 

where (W) is the average output power and (Q) the average heat flow from the QPC. 
One observes that a positive output power always contributes as a negative term to 


DQD-QPC: A stock, thermo, approach 


15 


the entropy production, which is compensated by the heat term (/3c — l^h) (Q) so as to 
satisfy the second law inequality (78). 

The efficiency p of the heat to work conversion process described above is defined 
by the ratio 




{W) -A/r(4) 


< Vc, 


(79) 


{Q) {Je) + {'^e) 

where rjc = l^h /fhe Carnot efficiency of the machine and the inequality is a direct 
consequence of the positivity of the entropy production (78). The Carnot efficiency may 
only be reached for a heat engine working reversibly, i.e. satisfying Si = 0. However, 
such machines work infinitely slowly so that the extracted power vanishes in the limit of 
a reversible machine. This issue has motivated the study of the efficiency at maximum 
output power in thermal engines [57, 58, 59, 60, 15, 61, 62]. In particular, the Curzon- 
Ahlborn efficiency 



has been shown to provide a universal upper bound on the efficiency at maximum power 
in machines working in the linear regime [58], as well as a good reference in the non¬ 
linear regime [61]. Reaching maximal output power and optimal conversion efficiency 
requires a fine tuning of the device parameters. For fixed A/r, the DQD spectrum 
and its coupling to the reservoirs via the tunneling amplitudes need to be adjusted. 
In particular, highest efficiencies are attained in the so-called regime of tight-coupling 
where the input and output power become proportional to each other [61, 62]. 

This can be understood in the present context by comparing the rates of the second 
order processes depicted in Fig. 4. In presence of the temperature gradient (77), the 
QPC will preferentially give energy to the DQD by inducing transitions from state |—) 
to state |-|-). The microscopic process leading to a net flow of charge against the bias 
A/i and involving such transitions is depicted in Fig. 4 a). Without optimization, this 
process is not more likely to happen than the other processes depicted in Fig. 4, which 
do not involve a charge transfer in the desired direction thus lowering the output power 
as well as the efficiency of the heat to work conversion. However, provided the tunneling 
rates satisfy 


bn-, F 2 -<C Fi_, r2+, (81) 

the process depicted in 4 a) becomes the dominant one, leading to a highly efficient 
conversion of the heat flow into electric output power. 

In the tight-coupling limit, i.e. for Fi_|_,r 2 - —?• 0, the matter current through the 
DQD and the energy current from the QPC are totally correlated, their average being 
thus proportional to each other 

{Je) + {Je) — 


(82) 
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Figure 4. Illustration of the different processes involving a transtion from state |—) 
to state 1+) in the DQD. 


as is verified by using the explicit expressions for the currents given in section 4.2. 
Consequently, the efficiency takes the simple form 


'/ — 

a;+- 

This last relation shows that in the tight coupling regime, our device only works as a heat 
engine producing a positive outpnt power in the range 0 < — A/r < a;+_. The similarity 
of expression (83) and the one obtained for the efficiency of the machine considered 
in Ref. [16] stems from the quantized character of the amount of energy exchanged 
between the reservoirs and the work converter in both models. In the present case, the 
QPC exchanges energy with the DQD in the form of quantas whose energy is given by 



UJ+-. 

In Fig. 5 a), the average output power is plotted against the rescaled efficiency 
p/pc for different values of the ratio 

. rv^iV 

Ti- r2+’ 

which measures the distance from the tight coupling regime (81). For each curve, the 
Carnot efficiency was held at the fixed valne pc = 0.6, whereas the DQD spectrum was 
numerically adjusted to maximize the output power at fixed A/r. Only values in the 
range where the device works as a heat engine are shown. As can be seen, the range 




















































Figure 5. a) Average output power as a function of the rescaled efficiency rijiric for 
different values of the asymmetry parameter d defined in (84). The Carnot efficiency 
was taken as 770 = 0.6. b) Efficiency at maximum power as a function of Carnot 
efficiency. The continuous red line is for Carnot efficiency while the dotted one follows 
the Curzon-Ahlborn efficiency (80). Remaining parameters were chosen for both plots 
as = I 5 ^0 = 0, Vs = V 4 = 0, (/ii+/i 2)/2 = 0.7, r 34 = r 43 = 0.01, Ti- = r 2 + = 0.1. 
The energies E+ and E- of the single-occupied states of the DQD, as well as the bias 
A/i in the case of the right plot, were numerically adjusted in order to reach a maximal 
output power. 


over which a positive output power can be produced as well as its magnitude decrease 
as one moves away from the tight coupling regime, i.e. as 0 increases. 

Fig. 5 b), shows curves of the efficiency where the DQD spectrum as well as the 
bias A/i were adjusted to reach the regime of maximum output power. As 9 increases, 
the efficiency is lowered due to the increasing contributions of the undesired processes 
described in Fig. 4. In the tight coupling regime, the thermodynamic efficiency attains 
values close to r]cA‘ 



Figure 6. a) Output power LDF. b) Stochastic efficiency LDF. Parameters are chosen 
as in Fig. 5 b) and the Carnot efficiency is set to rjc = 0.6 in both Figs. These PDFs 
thus characterize the output power and efficiency fluctuations at the intersection points 
between the dashed vertical line and the curves illustrated in Fig. 5. As explained 
in the text, the Carnot efficiency corresponds to the maximum of the LDF (dashed 
vertical line) and is thus the least likely to be observed in a single run experiment. 
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Up to this point we have considered the average properties of the output power 
and the so-called macroscopic efficiency, i.e. the efficiency defined as the ratio of the 
average output power over the average input power. In small devices displaying strong 
fluctuations, a more accurate characterization is provided by considering their statistical 
properties. Several works have recently shown that the so-called stochastic efficiency, 
defined along a single stochastic realisation of a thermal engine as 


>V 


(85) 



exhibits universal properties [41, 63, 64, 42, 65, 66]. We now briefly show that the PCS 
of the currents developed in Sec. 3 can be used in order to study these fluctuations. 

The probability distribution of a stochastic variable x = Ax/t is characterized 
in the long time limit by its large deviation function (LDP) X{x), that is, p{Ax) ~ 
exp{—for t ^ oo. It is shown in Ref. [42] that the LDPs for the output power 
and stochastic efficiency LDPs can respectively be obtained as 


X(>V) = max [Qw,q{Q^,0) — aW 
Xirjs) = -mmgy,^g{ar)s,a), 

a 


( 86 ) 

(87) 


a 


where the heat and work generating function gw,q{G(w,G(q) is obtained by setting 
^3 = ^4 = -^1 = —Apiayj and the other counting parameters to 0 in the current 

CGP (55). 

The LDPs of the fluctuating output power and efficiency are both illustrated on 
Pig. 6 a) and b) respectively, for different values of the asymmetry parameter chosen 
as in Pig. 5 and with the Carnot efficiency set to rjc = 0.6. Both output power and 
efficiency LDPs were evaluated in the regimes of maximal average output power. They 
thus characterize, respectively, the work fluctuations around the maxima’s of the curves 
depicted on Pig. 5 a), and the efficiency fluctuations at the intersection points between 
the curves and the dashed vertical line in Pig. 5 b). The minimum of the output 
power LDP, corresponding to the most probable output power in the long time limit, 
corresponds to the average value of the output power which increases as one gets closer 
to tight coupling, i.e. as 0 decreases. Similarly, the most likely value of efficiency in the 
long time limit, which lies the minimum of I{r)s), corresponds the macroscopic efficiency 
(79) evaluated in Pig. 5 b). Purthermore, one observes that its maximum (along the 
dashed vertical fine) lies at the Carnot efficiency {rjc = 0.6) in agreement with previous 
results [41, 63, 64, 42, 65]. The Carnot efficiency is thus the least likely to be observed 
in the long time limit. 

5.2. Regime of isothermal electric current conversion 

We now consider an isothermal regime (3 = [3j 'ij , in which the current in the electrically 
biased QPC is converted into work done against the electrical bias applied to the DQD. 
The DQD and QPC channels are subject to the chemical potential biases 


A/i = /ii — /i2 < 0 


0 < AjJiq = 10.3- ^^4, 


and 


( 88 ) 
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respectively. The irreversible entropy production in the system at steady state is then 
given by (see Eq. (74)) 

Si = ^(-{Wout)Pmn)) >0 (89) 

in terms of the input (Win) = ^Pq{JN) output (Wont) = power. 

The efficiency of the conversion process can be written in terms of the input and 
output powers as 

0<r]= < 1. (90) 

(y^in) 

A strong asymmetry in the tunneling amplitudes of the DQD channel was needed 
in order to achieve high power and efficiency in the regime of thermoelectric conversion. 
In the present case, a similar asymmetry in the tunneling amplitudes of the QPC is 
also mandatory as a consequence of the directional nature of the driving processes in 
the QPC. To understand this point, it is useful to consider the low temperature limit 
of the machine. In this limit, energy is transferred from the QPC to the DQD only 
if A/r^ > The several transfer processes involving a single electron transfer from 

reservoir 3 to 4 of the QPC, with associated transition in the DQD are depicted in Fig. 
7. By considering the rate constants involved in these two processes, we see that positive 
energy flow from the QPC to the DQD is enhanced provided 

P 43 “C Psr (91) 

so that the process depicted in Fig. 7 a) becomes the dominant one. 



Figure 7. Illustration of the processes involving a transfer of a single electron from 
reservoir 3 to 4 and the corresponding transition in the DQD system. 

On top of this, processes without net transfer of charge across the DQD like those 
depicted in Figs. 4 b) and c) are also undesirable since they waist energy. We thus 
assume both (91) and (81) to hold in order for our machine to work in optimal conditions. 

We further note that the DQD current and the contribution to the QPC current 
from electrons interacting with the DQD are tightly coupled in the regime where (81) 
and (91) are satisfied, that is (J^) oc {J%)d- However, due to the presence of bare 
tunneling events in the QPC, the DQD current and the total QPC current are not 
necessarily tightly coupled. Only in the case 7 (e) = 0 does the bare current in the QPC 
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vanish, {J%)d = 0, and the DQD current and the total current through the QPC become 
tightly coupled, i.e. (J^) oc {J%)d = (Jn)- 

This remark has important consequences on the properties of efficiency at maximum 
power, and illustrates well the crucial role played by the PCS formalism in determining 
all the thermodynamically relevant processes when analysing thermodynamic machines. 
The efficiency is illustrated in Fig. 8 for different values of the bare tunneling amplitude 
taken in the wide band limit, 7 (e) = 7 . We observe a significant difference in the 
qualitative behavior of efficiency depending on the value of the bare tunneling amplitude 
7 . In particular, we observe that 77 ^ 1 or 77 —?• 0 in the far from equilibrium regime, 
i.e. Afiq — 7 - 00 , depending on whether 7 = 0 or 7 7 ^ 0 respectively. 

The behavior of efficiency is best understood by writing it as 

—Apt, 

” ” (1 + S) ’ 

in terms of the bare current of electrons through the QPC {J%)i and where we assume 
both (81) and (91). 

For 7 = 0 , the bare current vanishes and the efficiency at maximum power in the 
equilibrium limit Afig ^ 0 converges to the value 77 = 1 / 2 , as predicted within linear 
response theory for systems working in the tight coupling regime [58]. In the large bias 
A fig — 7 - 00 , efficiency reaches the maximal value 77 = 1 . This is easily understood in the 
low temperature limit, i.e. /3 —>■ 00 . In this limit, the output power is non-zero provided 
that 



Afi < C(j+_ < Afdg. (93) 

When the bias Afig applied to the QPC is large, the values of Afi and a;+_ can be 
optimized to reach 
Apt, 


Aflg 


1. 


(94) 


For 7 7 ^ 0 , a fraction of the electrons tunneling through the QPC dissipates entropy 
without exchanging energy with the DQD thus lowering the efficiency as suggested by 
(92). Close to equilibrium, the efficiency converges to values below rj < 1 / 2 . 

To understand the properties of efficiency in the large bias limit Apig —>■ 00 we 
consider the behavior of the two contributions to the QPC current {Jff)d and in 

Fig. 9. We have assumed the tunneling amplitude 7 (e), to be a stepwise function equal 
to 7 on the interval [—a/2, o/2] and 0 elsewhere. In the wide band limit, i.e. o 00 , we 
see that the current (J|r)i is linearly growing and diverges as Ap,g ^ 00 as can be seen 
from its definition (65)-(66). This is in contrast to the current {Jff)d, which remains 
bounded VA/i^ since it is ultimately constrained by the DQD splitting as can be 
inferred from (65), (67) and the expressions for the rates (23) and (43). 

These remarks together with expression (92) for the efficiency and the fact that 
(J/f) remains bounded VA/ig explains why the efficiency decreases at least as fast as 
~ Api~^ in the large bias limit Afig ^ 00 as soon as 7 7 ^ 0 . 





DQD-QPC: A stock, thermo, approach 


21 



Figure 8. Efficiency of the current converter as a function of the input bias and for 
increasing bare tunneling parameter 7 . From top to bottom line we chose 7 = 10“^, 
5 X 10 -^ 10 -^ 5 X 10 -^ 10 -^ 5 X 10 -^, 10 “^, 5 x 10 “^ and 10 “^ The energies 
E+ and E_ and the bias A/i in the DQD were optimized to reach maximum output 
power. The remaining parameters were chosen as = 1 , (El + E 2)/2 = 1 . 5 , Eq = 0 , 

r2+ = ri_ = r34 = o.i. 


a) b) 






— a=oo 
T a=l 
^ a=2 
-f- a=3 


Figure 9. Curves of the current contributions in the QPC which do (left plot) and do 
not (right plot) induce transitions in the DQD, for different values of the band width 
a (see text). Parameters have been chosen as /3 = 1, Ei = E 2 = 1-5, Eq = 0, E+ = 0.5, 

E- = 0.1, r 2 + = Pi- = r34 = 7 = 0 . 1 , 


We have thus shown that the bare tunneling events in the QPC can highly reduce 
the machine efficiency. The PCS formalism is here crucial in order to keep track 
of thermodynamically relevant processes which are otherwise missing in a stochastic 
description of the DQD populations. Let us finally mention that such processes did not 
affect the performance of the thermal engine considered in the present section due to the 
fact that the matter current in the QPC does not contribute to the entropy production 
since its corresponding affinity was set to zero, i.e. /is = /i 4 . 

6. Conclusion 

We fully characterized the nonequilibrium thermodynamics of a circuit composed of a 
double quantum dot (DQD) and quantum point contact (QPC) channels within the 
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framework of stochastic thermodynamics. 

By using the modified quantum master equation formalism, we identified and 
provided a detailed description of all the microscopic processes contributing to the 
entropy changes in the system and the environment. We showed that the transition rates 
of processes related by time reversal satisfy the local detailed balance condition. This 
condition holds for the rates of electron transfers between the DQD and its reservoirs 
as well as for the rates of energy exchange processes with the ont-of-eqnilibrinm QPC. 
We also established a steady-state finctnation theorem (FT) for the heat and matter 
currents across both channels, which reduces to a bivariate FT for the matter currents 
across each channel in an isothermal circuit. 

Based on this analysis, we considered two regimes where the circuit operates as 
a thermoelectric or electric converter. We identified the optimal working condition in 
both cases and evaluated the statistics of output power and efficiency at steady state. 

Our study illustrates very well how stochastic thermodynamics enables one to 
structure the analysis of the transport properties of non-trivial mesoscopic devices such 
as the DQD-QPC circuit and to discriminate the universal (i.e. thermodynamic) features 
from the system specific ones. This theory has become an essential tool for quantum 
transport in the weak coupling regime. 
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Diagonalisation of the DQD Hamiltonian 


Diagonalisation of the DQD Hamiltonian (3) can be performed analytically. We first 
introduce the basis of single dot occupation states as 


IUOb) = c^|0^0b) 

(95) 

|0 a1b) = ctlOAOfi) 

(96) 

where 0^0^) denotes the empty state of the DQD. 

The eigenstates of the DQD 

Hamiltonian (3) are then given by 

|0 ) = IO^Ob) 

(97) 

1 +) = cos^IUOb)+ sin^|0^1s) 

(98) 

9 9 

|-) = sin-|lA0s) -cos-|0^1s) 

(99) 

in terms of the mixing angle defined by 

^ 2r 

tan 9 = -. 

(100) 
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The corresponding energies are given in terms of the localised basis parameters by 

Eo = 0 E±= ± (101) 

Other parameters of the Hamiltonian in the single and manybody basis are related by 
T^+ = t^A cos ^ (102) 

= (103) 

^2+= 4b sin ^ (104) 

^2- = - 4b cos (105) 

and 

rS = l(tf + tf ) ±\(tf - 4‘')cos9 (106) 

ytt' ^ yU' ^ K^kk' _ 



